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^ ■ Abstract 
T— I ■ 

^ ' There are two processes of matter creation after inflation that may be relevant to 

CN . the resolution of the puzzle of cosmic rays observed with energies beyond GZK 
cut-off: 1) gravitational creation of superheavy (quasi)stable particles, and 2) non- 

^f) ' thermal phase transitions leading to formation of topological defects. We review 

O ■ both possibilities. 

r~| ' Key words: Cosmic rays; Dark matter; Early universe 

Ph: 

> 
X 

• Cosmological and astrophysical considerations are able to provide the strongest 

restrictions on parameters of particle physics models and even rule out some 
classes of models entirely. This is especially valuable when the model is unre- 
stricted by laboratory experiments (which is often the case). Among famous 
results which made a strong impact on model building is the cosmological 
domain wall problem which appears in models with spontaneous breaking of 
discrete symmetries [1] and the problem of magnetic monopoles in Grand Uni- 
fied theories [2]. In return, studies of cosmological phase transitions [3] and 
of the dynamics of bubbles of a metastable vacuum [4] lead to the change 
of basic concepts of the cosmology of the early Universe, and inflationary 
cosmology [5-8] was born (for reviews see [9,10]). Inflation gives a possible 
solution to horizon, flatness and homogeneity problems of "classical" cosmol- 
ogy [6]. Inflation was designed to solve the problem of unwanted relics, like 
magnetic monopoles. It was promptly realized [11] that inflation can generate 
small amplitude large scale density fluctuations which are the necessary seeds 
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for the galaxy and the large scale structure formation in the Universe. This 
elevates inflation from the rank of a "broad brush problem solver" into the 
rank of a testable hypotheses. And testable in fine details, as rapidly accu- 
mulating data on cosmic microwave background fluctuations (starting from 
COBE detection [12] through numerous balloon and ground based CMBR ex- 
periments and with culmination at MAP [13] and PLANK [14] anticipated 
detailed maps of anisotropy of the microwave sky) and huge galaxy catalogs 
like the already collecting data SLOAN digital sky surview [15] will provide a 
wealth of cosmological information. 

Inflation is generally assumed to be driven by the special scalar field known 
as the inflaton. During infiation, the infiaton field slowly rolls down towards the 
minimum of its potential. Inflation ends when the potential energy associated 
with the inflaton fleld becomes smaller than the kinetic energy, which happens 
when magnitude of the inflaton fleld decreases below the Plank scale, < 
Mpi and "cold" coherent oscillations of the inflaton field commence. These 
oscillations contained all the energy of the Universe at that time. 

All matter in the Universe was created in reheating, which is nothing but 
decay of the zero momentum mode of inflation oscillations. The process is 
obviously of such vital importance that here too one may hope to flnd some 
observable consequences, specific for the process itself and for particular mod- 
els of particle physics, despite the fact that scales relevant for the reheating 
are very small. And, indeed, we now believe that there can be some clues left. 
Among those are: topological defects production in non-thermal phase transi- 
tions [16], GUT scale baryogenesis [17], generation of primordial background 
of stochastic gravitational waves at high frequencies [18], just to mention a 
few. However, the most interesting can be a possible relation to a mounting 
puzzle of the Ultra High Energy Cosmic Rays (UHECR) [19]. 

When a proton (or neutron) propagates in CMB, it gradually looses energy 
colliding with photons and creating pions. There is a threshold energy for 
the process, so it is effective for very energetic nucleons only, which leads to 
the Greisen-Zatsepin-Kuzmin (GZK) cutoff [20] of the high energy tail of the 
spectrum of cosmic rays. All this means that detection of, say, 3 x 10^° eV 
proton would require its source to be within ~ 50 Mpc. However, several well 
established events above the cut-off were observed by Yakutsk [21], Haverah 
Park [22], Fly Eye [23] and AGASA [24] collaborations (for the recent reviews 
see Refs. [25,26]). 

Results from the AGASA experiment [27] are shown in Fig. 1. The dashed 

curve represents the expected spectrum if conventional extragalactic sources 
of UHECR would be distributed uniformly in the Universe. This curve displays 
the theoretical GZK cut-off, but one observes events which are way above it. 
(Numbers attached to the data points show the number of events observed in 
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Fig. 1. AGASA data set [27], February 1990 - October 1997. 

each energy bin.) Note that no candidate astrophysical source, hke powerful 
active galaxy nuclei, were found in the directions of all six events with E > 10^^ 
eV [27] (at these energies cosmic rays experience httle deflection by galactic 
magnetic fields). 

Is some unexpected astrophysics at work here or this is at last an indication 
of the long awaited new physics ? 

There are two logical possibilities to produce UHE cosmic rays: either charged 
particles have to be accelerated to energies E > 10^° eV, or UHECR originate 
in decays of heavy X-particles, nix > ^0^^ GeV. Maximum energy which can 
be achieved in an accelerating site of the size R which has the magnetic field 
strength B is [28]: 

E<10''z4^J-eV . (1) 
fxG Mpc ^ ^ 



A magnetic field is required either to keep the particle confined within the ac- 
celerating region or to produce an accelerating electric field. For protons (Z=l) 
a few sources satisfy this condition: pulsars, active galactic nuclei (AGNs) and 
radio-galaxies. However, energy losses (pair production and meson photopro- 
duction) restrict the maxim energy to ii^ < 10^^ eV in pulsars and AGNs 
[25,29], while radio-galaxies that he along the arrival directions of UHECR 
are situated at large cosmological distances, > 100 Mpc [30], i.e. beyond the 
GZK radius. Similar conclusion seem to be true with respect to cosmological 
Gamma Ray Bursts as a possible source of UHECRs [31]. 

New astrophysics which may work is a possibility to generate UHECR within 
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GZK sphere in remnants of dead quasars [32] (these are dormant galaxies 
which harbour supermassive spinning black hole). 

New physics suggested as an explanation of UHE cosmic rays, range up to the 
violation of the Lorentz invariance [33] . Among less radical extensions of the 
standard model are: 

• The existence of a particle which is immune to CMB in comparison with 
nucleons. In this scenario the primary particle is produced in remote astro- 
physical accelerators (e.g. radio-galaxies) and is able to travel larger cos- 
mological distances while having energies above the GZK cut-off. There are 

variations to this scheme. 

• Suprsymmetric partner of gluon, the gluino, can form bound states with 
quarks and gluons. If gluino is hght and quasistable (see e.g. [34,35]), 
the hghtest gluino containing baryons will have sufficiently large GZK 
threshold to be such a messenger [34] and as a hadron it will be able to 
produce normal air showers in the Earth's atmosphere. However, there are 
strong arguments due to Voloshin and Okun [36], against light quasistable 
gluino based on constraints on the abundance of anomalous heavy isotopes 
which also will be formed as bound states with gluino. 

• High energy (anti)neutrinos produced in distant astrophysical sources will 
annihilate via Z'^ resonance on the relic neutrinos and produce energetic 
gamma or nucleon [37]. The relic neutrino masses in the eV range are 
consistent with this scenario [38], as well as with the Super-Kamiokande 
results. The required high density of the relic neutrinos is achieved if gravi- 
tational clumping takes place [37] or if the Universe has a significant Icpton 
asymmetry in background neutrinos [38]. Even then the total luminosity 
of the neutrino sources in the Universe must be as high as 10^^ — 10^ of 
its photon luminosity, and, therefore, neutrino-only sources are called for 
by the upper bound from the flux of the cosmic rays [39]. An indepen- 
dent constraint on the density of the relic neutrinos comes from CMBR 
and already the present data start to be challenging for models with large 
neutrino asymmetry [40]. 

• Another class of suggestions is related to topological defects. UHECR are 
produced when topological defects decompose to constituent fields (X-particles) 
which in turn decay [41]. Maximum energy is not a problem here, but in 
models which involve string [42] or superconducting string [41] networks, 
the typical separation between defects is of order of the Hubble distance 
and thus these models are subject to GZK cut-off. Models in which defects 
can decay "locally" include networks of monopoles connected by strings 
(necklaces) [43], vortons (charge and current carrying loops of supercon- 
ducting strings stabilized by angular momentum) [44], and monopolonium 
(bound monopole-antimonopole pairs) [45]. Finally, magnetic monopoles ac- 
celerated by intergalactic magnetic fields were also considered as primary 
UHECR particls [46]. 
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Fig. 2. Predicted fluxes from decaying X-particles, as calculated in Ref. [49] and the 
data. Latest AGASA results, Fig. 1, are not shown. 

• Conceptually the simplest possibility is that UHECR are produced cos- 

mologically locally in decays of some new particle [47,48]. GZK cut-off is 
automatically avoided but the candidate X-particle must obviously obey 
constraints on mass, number density and lifetime. 



2 UHECR from decaying particles 

In order to produce cosmic rays in the energy range E > 10^^ GeV, the 
decaying primary particle has to be heavy, with the mass well above GZK cut- 
off, mx > 10^^ GeV. The lifetime, rx, cannot be much smaller than the age of 
the Universe, % 10^° yr. Given this shortest possible hfetime, the observed 
flux of UHE cosmic rays will be generated with the rather low number density 
of X-particles, Qx ~ 10~^^, where Qx = ^x^x/ Pent-, ^x is the number 
density of X-particles and Pcrit is the critical density. On the other hand, 
X-particles must not overdose the Universe, ilx < 1- With Qx ~ 1, the X- 
particles may play the role of cold dark matter and the observed flux of UHE 
cosmic rays can be matched if tx ~ 10^^ yr. 

Spectra of UHE cosmic rays arising in decays of relic X-particles were success- 
fully fitted to the data for mx in the range 10^^ < mx/GeV < 10^^ [49,50]. 
For example, the fit of Berezinskii et. al. [49] to observed fiuxes of UHECR as- 
suming mx ~ 10^^ GeV is shown in Fig. 2. Beside the mass of the X-particle 
there is another parameter which controls the flux of the cosmic rays from 
decaying particles: namely, the ratio of X-particles number density and their 
lifetime. For the fit in Fig. 2 it was used {^x /^CDM){tu /tx) — 5 x 10~^^. 
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The problem of the particle physics mechanism responsible for a long but 
finite lifetime of very heavy particles can be solved in several ways. For exam- 
ple, some otherwise conserved quantum number carried by X-particles may 
be broken very weakly due to instanton transitions [47], or quantum gravity 
(wormhole) effects [48]. If instantons are responsible for X-particle decays, 
the lifetime is estimated as tx ~ m'^ ■ exp(47r/Q;x), where ax is the coupling 
constant of the relevant gauge interaction. The lifetime will fit the allowed 
window if the coupling constant (at the scale mx) is ax ~ 0.1 [47]. 

A class of natural candidates for superheavy long-living particles which arise in 
string and M theory was re-evaluated recently in Refs. [51] and particles with 
desired mass and long life-time were identified. Other interesting candidates 
were found among adjoint messengers in gauge mediated supergravity models 
[52] and in models of superheavy dark matter with discrete gauge symme- 
tries [53] . Superheavy dark matter candidates in superstrings and supregarvity 
models were considered also in Refs. [54]. 

Below we address the issue of X-particle abundance. 



3 Superheavy particle genesis in the eeirly Universe 

Seperheavy particles can be created in the early Universe by several mecha- 
nisms. Among those are: 

• Non-equilibrium "thermal" production in scattering or decay processes in 
primordial plasma [47,48]. 

• Production during decay of inflaton oscillations ("preheating") [55-58]. 

• Direct gravitational production from vacuum fluctuations during inflation 
[59,19,60]. 

In any case the final ratio of the density in X particles to the entropy den- 
sity is normalized by the reheating temperature. The reheating temperature 
is limited to the value below 10^ - 10^ GeV in supergravity models with de- 
caying heavy gravitino [61]. This restricts model parameters when "thermal" 
mechanism of heavy particle production is operative (but does not rule it out 
[47,48,62]). 

The last two mechanisms are closely related to each other and both can be 
described on equal footing within frameworks of a single unified approach: 
particle creation in external time varying background. However, while the 
outcome of the second mechanism is highly dependent upon the strength of 
the interaction of the X-field to the infiaton, no coupling (e.g. to the infiaton 
or plasma) is needed in the third mechanism when the temporal change of 
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the metric is the single cause of particle production. Even absolutely sterile 
particles are produced by the third mechanism which may be relevant for very 
long-living superheavy particles. Resulting abundance is quite independent of 
the detailed nature of the particle which makes the superheavy (quasi)stable 
X-particle a very interesting dark matter candidate. 

We concentrate here on the second and third mechanisms and from the start 
we introduce coupling of the X- field to the inflaton for uniformity of discussion. 
The limit of zero coupling will correspond to pure gravitational production. 

In the case of a heavy scalar field X we consider the model 

L = \{d,4>f - v{<t>) + \{d,xf - \{Ml - CR)X' - ^0^X2 . (2) 

Here V{(f)) is the inflaton potential. In simple "chaotic" [8] models of inflation 
the inflaton is either "massive" with the scalar potential V{(f)) = m^0^/2, or 
"massless", V{(p) = A0*^/4. Normalization to the large scale structure requires 
m|/Mpi 10~^^ in the former model and A ^ 10~^^ in the latter model. 
The constant ^ describes direct coupling to the space-time curvature R, with 
^ = corresponding to the minimal coupling and ^ = 1/6 being the case of 
conformal coupling. 

Fermion field (spin |) is conformally coupled to gravity. In addition to stan- 
dard kinetic and mass terms it also may have coupling to the infiaton, Vi — 

gct>XX. 

It is convenient to work in conformal metric ds^ = a{ri)'^{dTf — rfx^) with 
rescaled fields, tp = 4>a{rj) and % = Xa{rjY , where s = 1 and s = | for scalar 
and fermion fields respectively. In what follows we measure time and space 
intervals in units of infiaton mass, r = mr]. 

3.1 Quantum fields in classical backgrounds 

Here we summarise the basic formalism of particle creation in external classical 
background (e.g. space-time metric of an expanding universe or oscillating 
infiaton field). For more details see e.g. Refs. [63-66]. 

i) Spin bosons. 

A real scalar field is Fourier expanded in a comoving box 

X{r, x) = Y.\Xk{r)a^. + xl{r)al^¥'^- . (3) 

k 



7 



Annihilation and creation operators commute except for [ak-oy = 1. The 
mode functions, Xk = Xki^) of ^ scalar Bose field are solutions of the oscillator 
equation 

x\ + col{r)xk^O: (4) 



with the time dependent frequency 

ulir) ^k' + -{6^-l) + my + Aq^' , (5) 



a 



where ' = d/dr and 



Arm? 



Here 0(0) is the value of the inflaton field when it starts to oscillate (which 
corresponds to normalization (^(0) = 1). 

Let ujk at some time interval satisfy the adiabatic condition <^ 1. 

We can choose solutions of Eq. (4) which enter decomposition Eq. (3) to be 
positive- frequency modes Xk — ^k e~*'^'='^ and define vacuum state ak|0) = 0. 
The number of particles in a non-vacuum state will be constant during evo- 
lution through this adiabatic interval. Let the adiabatic condition be violated 
for some time and then the system enters another adiabatic interval. In that 
interval we can define another set of positive- frequency modes and correspond- 
ing vacuum. Initial positive-frequency modes evolved through non-adiabatic 
region will not coincide with "out" state modes, but one set of modes can 
be expressed in terms of the other. This decomposition is called Bogolyubov 
transformation [67] . Since one and the same field is expanded with the use of 
two different sets of mode functions, the Fourier coefficients are also related 
to each other 

al'' = atak + /5^aj, . (7) 

It follows immediately that the initial vacuum state at late times contains 
particles 

{Q\al°-'ar\0) = W . (8) 



Technically it is easier to find Bogolyubov's coefficients by diagonalizing Hamil- 
tonian of the field X. For any time moment r this procedure gives 

|^j, lx-.P+^^y-2. _ (9) 
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where mode functions are solutions of Eq. (4) with initial (vacuum) conditions 

Xk{0)=u^'/', x'k{0) = -tioXk- (10) 



ii) Spin 1/2 fermions. 

The relevant mode functions of the Fermi field satisfy the oscillator equation 
with the complex frequency 

x"k + {^l-^KE)xk = o , (11) 



where the real part of the frequency is given by ul = k"^ + mlg and mefr = 
m^a + ^/q(p. We choose 



X,(0) = Jl - x'M^-i^Xk, (12) 

V cu 



as the initial conditions. In this case we find per spin state 



2 uj-rries- Im{xkXk) 



l^'l- -J—'' . (13) 



Finally, the number density of X-particles created by time varying background 
is 



1 

nx 



J2 J m'k'dk , (14) 



where J2s is the sum over spin states. The expression (14) gives the number 
density of particles only, with an equal amount of antiparticles being created 
in the case of charged fields. 



3.2 Gravitational creation of particles 



It was noticed [59,19] that superheavy particles arc produced gravitationally 
in the early Universe from vacuum fluctuations and their abundance can be 
correct naturally, if the standard Friedmann epoch in the Universe evolution 
was preceded by the inflationary stage. This is a fundamental process of par- 
ticle creation unavoidable in the time varying background and it requires no 
interactions. Temporal change of the metric is the single cause of particle pro- 
duction. Basically, it is the same process which during inflation had generated 
primordial large scale density perturbations. No coupling (e.g. to the inflaton 
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or plasma) is needed. All one needs are stable (very long-living) X-particles 
with mass of order of the inflaton mass, rrix ~ 10^^ GeV. Inflationary stage 
is not required to produce superheavy particles from the vacuum. Rather, 
the inflation provides a cut off in excessive gravitational production of heavy 
particles which would happen in the Priedmann Universe if it would start 
from the initial singularity [19]. Resulting abundance is quite independent of 
the detailed nature of the particle which makes the superheavy (quasi) stable 
X-particle a very interesting dark matter candidate. So, we start our consider- 
ation with gravitational creation of particles, i.e. we put g = {or equivalently 
q — 0) in the formulas above. 

3.2.1 Friedmann Cosmology 

For particles with conformal coupling to gravity (fermions or scalars with 
^ = 1/6), it is the particle mass which couples the system to the background 
expansion and serves as the source of particle creation. Therefore, just on 
dimensional grounds, we expect 

nx oc mx(i~^ (15) 

at late times when particle creation diminishes. In Friedmann cosmology, a oc 
(mt)" oc {m/H)" {a — j and o; = | for radiation and matter dominated 
expansion respectively). We conclude that the anticipated formulae for the 
X-particles abundance can be parameterised as 



On the other hand, it is expansion of the Universe which is responsible for 
particle creation. Therefore, this equation which describes simple dilution of 
already created particles, is valid when expansion becomes negligible, H <C 
mx- This means also that particles with nix ^ H cannot be created by this 
mechanism. Creation occurs when H ~ mx- The coefficient Cq, depends upon 
the background cosmology only, and it can be found numerically [19], see 



In particular, for the radiation dominated Universe, which was studied also 
in Ref. [63], one finds Qx = Px/Pc = mxnx ^'^'^Gt'^ /3 with the present value 
of Qx being equal to Qx ~ 2 x 10~'^{mx/Mp^)^/rnxt^^ where tg is the time 
of equal densities of radiation and matter in the n — 1 Universe. This gives 
rix ~ (mx/lO^GeV)^/^. Stable weakly interacting particles with mx ~ 10^ 
GeV will overdose the Universe even if initially they were in a vacuum state 
and were created from the vacuum during the regular radiation dominated 
stage of the Universe evolution. 




(16) 



Fig. 3. 
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Fig. 3. The coefficient Cq, defined in Eq. (16), is shown as a function of a for the 
background cosmology with a power law scale factor a oat"', Ref . [60] . 

There is no room for superheavy particles in our Universe if it started from 
the initial Priedmann singularity [19], since the value of the Hubble constant 
is limited from above only by the Planck constant in this case. 



3.2.2 Inflationary Cosmology 

However, this restriction will not be valid if inflation separates initial con- 
ditions, whatever they were, from the observable Universe. In inflationary 
cosmology the Hubble constant (in effect) did not exceeded the inflaton mass, 
H < m^. The mass of the inflaton field has to be ~ 10^^ GeV as con- 
strained by the amplitude of primordial density fluctuations relevant for the 
large scale structure formation. Therefore, direct gravitational production of 
particles with mx > ~ 10^^ GeV has to be suppressed in inflationary 
cosmology. 

Particle creation from vacuum fluctuations during inflation (or in the de Sitter 
space) was extensively studied [68,69], usually in the case of small mx and in 
application to generation of density fluctuations necessary for the large scale 
structure formation. The characteristic quantity which is usually cited in this 
applications, the variance of the field, {X"^) , is defined by an expression similar 
to Eq. (14), in the typical case ctfc ~ —13k the integrand is being multiplied 
by the factor 2siv?{uJhT)/uJk. Por example, for the scalar Bose field with the 
minimal coupling to the curvature, (X^) = SHf'/Sn'^m'j^ if mx ^ i?j[68)69]. 
Por massless self-interacting field (X^) ^ 0.132H^ /\/\ [70]. Particle creation 
for the specific case of the Hubble dependent effective mass, mx{t) oc i^(t), 
was considered in Ref. [71]. 
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Fig. 4. Ratio of the energy density in X-particles, gravitationally generated in infla- 
tionary cosmology, to the critical energy density is shown as a function of X-particle 
mass, Ref. [60]. 

Results [60] of direct numerical integration of gravitational creation of super- 
heavy particles in chaotic inflation model with the potential = m00^/2 
is shown in Fig. 2. 

This figure was calculated assuming Tr = 10^ GeV for the reheating temper- 
ature. (At reheating the entropy of the Universe was created in addition to 
X-particles. In general, multiply this figure by the ratio GeV and devide 

it by the fractional entropy increase per comoving volume if it was significant 
at some late epoch.) The reheating temperature is constrained, Tr < 10^ GeV, 
in supergravity theory [61]. We find that Qxh'^ < 1 if mx ~ (few) x 10^^ GeV. 
This value of mass is in the range suitable for the explanation of UHECR 
events [19]. Gravitationally created superheavy X-particles can even be the 
dominating form of matter in the Universe today if X-particles are in this 
mass range [59,19]. 



3.2.3 Isocurvature fluctuations in superheavy particle matter 

In numerical calculations, Ref. [60], it was found that the variance (X^) of the 
field X measured at the end of infiation is independent upon mx if the mass 
of X is small and coupling to curvature is minimal. At some later epoch when 
H « mx (which will be long after the end of inflation if X is a fight field) 
the field X starts to oscillate on all scales, including k = 0. Only at this time, 
which we denote by tx, all field fiuctuations are transformed into the non-zero 
particle density and we can use px = thxtix ~ The variance of X 

fluctuations was unchanged on large scales, starting from the end of inflation 
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Fig. 5. Spectrum of created particles, k^n{k), in a model with massive inflaton is 
shown for several choices of the mass of scalar X-particle with the minimal cou- 
pling (solid lines) and the conformal coupling (dotted line), Ref. [60]. Masses and 
momenta, k, are given in units of the inflaton mass. 

down to the time tx- So, when the field starts to oscillate px oc m|-. However, 
the energy density of the inflaton field, p — SH^/SnG, decreased during this 
time interval in proportion to H'^{tx) / H'^{0) ~ m\/ H'^iQ). That is why the 
ratio of the energy density in X-particles to the total energy density does not 
depend on mx when measured at t > tx, see Fig. 4. 

Variance of the field X is different from the usually calculated for the fixed 
de Sitter inflationary background because we consider the actual evolution 
of the scale factor and the value of the Hubble parameter is not constant 
during inflation, being larger at earlier times. Correspondingly, the number 
of created particles per decade of k grows logarithmically towards small k if 
mx is small. (Power spectrum behaves similarly.) The examples of the particle 
number, k^nx{k). for several values of mx are shown in Fig. 5 at the moment 
corresponding to 10 completed inflaton oscillations. The particle momentum 
is measured in units of the inflaton mass. In contrast to this, in the flxed de 
Sitter background A7^'^{X'^) ^ j' d\nk [k / Hf-^^ with (3-2z/) ^ 2m\l?>H'^ 
at small mx, and consequently {X"^) oc l/m\. Note that the power spectrum 
in the flxed de Sitter background grows towards large values of A;, which is 
opposite to the behaviour of Fig. 5. 

Therefore, calculations which would be based on the customary procedure of 
matching a fixed de Sitter background to a subsequent Friedmann stage would 
give wrong results, with ^ oo at mx — 0. 

Matching is also dangerous in the case of large mx- When the change is too 
abrupt, it generates artificial particles. This may easily happen for mx > iri^, 
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see e.g. [72] where excessive production was found. At mx ~ the number 
of created particles decreases exponentially with nix- 

As Fig. 5 shows, the power spectrum of fluctuations in X-particles is almost 
scale independent at small k if mx/^^ ~ 1- Therefore, if such particles 
constitute a considerable fraction of dark matter, these fluctuations will be 
transformed into isocurvature density perturbations at late times and can af- 
fect large scale structure formation. Isocurvature fluctuations produce 6 times 
larger angular temperature fluctuations in cosmic microwave background ra- 
diation (CMBR) for the same amplitude of long-wavelength density pertur- 
bations compared to the adiabatic case [73]. Therefore, to fit observations by 
a single spectrum, the mass fluctuation spectrum in isocurvature cold dark 
matter cosmology must be tilted (with respect to scale invariant spectrum) to 
favor smaller scales. 

Let the power spectrum of the field fluctuations be k^Px{k) oc k^. Fit to the 
second moments of the large-scale mass and cosmic microwave background 
distributions requires f3 > 0.25. Models with (3 ranging from 0.3 to 0.6 were 
considered in Refs. [74,75]. It is interesting that the spectrum is indeed cor- 
rectly tilted for mx which is somewhat larger than the inflaton mass, see 
Fig. 5. On the other hand we see that lighter particles, mx < rn^, with mini- 
mal coupling to gravity and VLx ~ 1 are excluded. 



3.3 Preheating 

Considering gravitational creation of particles only, we would be left today 
with an oscillating inflaton field dominating the Universe.^ To reheat the 
Universe we must couple the inflaton to some other fields. In oscillating back- 
ground and in flat space-time the number density of particles grows exponen- 
tially: parametric resonance is always effective, see e.g. Ref. [66]. This would 
lead to explosive decay of oscillations. However, in an expanding universe the 
resonance is blocked by the redshift, which removes created particles out of 
the resonance bands. Thus, coherent oscillations of e.g. axion [77], or moduli 
fields [78], do not decay via parametric resonance in the expanding Universe. 
It is interesting that the field trapped in a (self-) gravitational well can burst 
in radiation, in principle [79,80], but limitations on the relation between den- 
sity and size of the clump make it hard to achieve critical conditions when 
elaborated dark matter models, like the axion [79,81] are considered (axion 
miniclusters are promising objects in this respect though [82]). 



Exclusions from this rule exist if inflaton potential does not have a minimum and 
energy density in the inflaton fleld after inflation decreases faster than in conven- 
tional matter [76]. 
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Narrow width of the resonance band is a prime reason of stabihty and such 
fields whill decay only when the Hubble constant H becomes smaller than the 
particle width F (i.e. the life-time of individual particle becomes smaller than 
the age of the Universe). This results in reheating temperature Trh ~ \/TMpi. 
For a long time it was believed [83] that this is the end of the story for the 
inflaton field, until it was realized [55] that the resonance parameter g, Eq. (6) 
can naturally be extremely large in the inflaton caseQ. Indeed, the coupling 
constant is multiplied by the ratio of initial amplitude of the field to its mass, 
squared, which for the inflation is ~ 10^°. The resonance is broad and it is 
impossible to red-shift the system out of the resonance. Resulting explosive 
decay is not actually a parametric resonance, since resonance parameter can 
change by orders of magnitude just during one period of oscillation [56,85], 
but decay occurs anyway because each period of oscillations the adiabaticity 
conditions are violated during some short time intervals [56,85]. It is possible 
that the inflaton oscillations decay after only a dozen of oscillations. 

Let us now consider this process and introduce a non-zero coupling of the 
inflaton and the X fields. For the sufficiently large value of q in Eq. (6) the 
gravitational production of particles is negligible. If the amount of produced 
particles is relatively small, the process of their creation still can be considered 
as creation by external time varying background (oscillating zero mode of 
the inflaton). However, production can be very efficient at large q and back 
reaction of produced particles on the motion of the inflaton field has to be 
included. This restores conservation laws and, as a matter of fact, we are 
considering the decay of the inflaton into X now. Here we restrict ourselves 
to the case of scalar X-particles. 

The equation of motion for the inflaton fleld is 



Now the background (i.e. the inflaton fleld) is not homogeneous anymore and 
the complete quantum problem is not easy. However, the efficiency of parti- 
cle creation saves the day: the system rapidly becomes classical [86] and the 
problem can be handed over to a computer [56]- [58]. 

Instead of the number density of created particles it is more convenient (and 
rigorous) to measure the related quantity: the variance of the fluctuation fleld 



For a discussion of the parametric resonance induced by an inflaton field see also 



ip" - VV ^ + + 4gxV = • 



(17) 



a 




(18) 



refs. [84]. 
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Initially fluctuations are small and the problem can be followed in the Hartrcc 
approximation. Namely, the field x in Eq. (17) is replaced by its average, (x^), 
and spatial gradients of the inflaton field are neglected. The formalism of sec- 
tion 3.1 then applies. In the situation of efficient particle creation occupation 
numbers grow and at some point quantum averages can be approximated by 
classical averages computed with the help of a certain distribution function 
[87,86]. At late times the problem becomes classical and it has to be supple- 
mented with appropriate initial conditions which refiect early quantum evo- 
lution. Those are specified in the Fourier space as follows [86,56]. Amplitudes 
of mode functions are distributed with the probability density 



P[Xk\ ^ exp 



■m.(o)ix.(o)r' 



(19) 



and mode functions have random phases. The initial "velocities" are locked 
to their corresponding "coordinates" 

X,-(0) = -iu^mXk ■ (20) 



The parameter m^/0^(O) sets the scale of (x^) which separates regions of 
quantum and classical fiuctuations in this model. The semiclassical description 
is rehable as long as g'^/47r = qm'^/A7r(f)'^{0) <^ 1. 

Full non-linear problem can now be solved numerically on the lattice (namely, 
Eq. (17) and corresponding equation for the field x(r, x) are solved directly 
in the configuration space). This classical problem accounts for all the ef- 
fects of particle creation, their rescattering, inverse decays, etc. Initial stage 
of intensive growth of fluctuations (similar to parametric resonance, but not 
equivalent to it since the resonance parameter q changes rapidly in the ex- 
panding Universe) is followed by the stage which is similar to the Kolmogorov 
turbulence when smooth power spectra of fluctuations are established which 
slowly approach equilibrium [86,57,58,88]. 

An important quantity is the maximum strength of fluctuations achieved dur- 
ing the course of evolution (fluctuations are diluted by the expansion of the 
Universe after decay slows down). Compilation of results of Refs. [56,57] for 
the realistic case m^;, = lO^^Afpi is shown in Fig. 6 as a function of model pa- 
rameters. The stars are results of the lattice calculation which takes all back 
reaction effects into account, solid lines correspond to the computationally less 
expensive Hartree approximation. 

To compare this with the gravitational creation note that in the Hartree ap- 
proximation px/p ~ 2 when fluctuations reach their maximum. Even the 
lowest level of fluctuations shown in Fig. 6 exceeds the amount of gravitation- 
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Fig. 6. The maximum value of the variance (^max) ^ ^ function of q. Sohd hues 

correspond to the Hartree approximation with different values of = nix/nT-tf,, 
Ref. [56]. The stars are results of fully non- linear lattice calculations with mx ^ m^, 
Ref. [57]. 

ally created particles. However, while stable particles will be generally over- 
produced, this mechanism still can be relevant for dark matter and UHECR 
phenomenon, when e.g. for a given value of nix the value of q is low enough to 
prevent effective particle creation, but not negligible. Fig. 6 shows that parti- 
cles 10 times heavier than inflaton are produced and even this is not a limit 
in the case of very large q [89]. 

In the opposite situation when q is sufficiently large and fluctuations develop 
to the full strength, an interesting and important phenomenon of non-thermal 
phase transitions [16] can occur which is the subject of discussion in the fol- 
lowing section. 



4 Topological defects and inflation 

Decaying topological defect can naturally produce very energetic particles, 
and this may be related to UHECR [41-46], for recent reviews see [26]. How- 
ever, among motivations for inflation there was the necessity to get rid of 
unwanted topological defects. And inflation is doing this job excellently. Since 
temperature after reheating is constrained, especially severely in supergravity 
models, it might be that the Universe was never reheated up to the point 
of GUT phase transitions. Topological defects with a sufficiently high scale 
of symmetry breaking cannot be created. Then, how could such topological 
defects populate the Universe? 
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The answer may be provided by non-thermal phase transitions [16] which 
can occur in preheating after inflation. Explosive particle production caused 
by stimulated decay of inflaton oscillations lead to anomalously high field 
variances which restore symmetries of the theory even if the actual reheating 
temperature is small. The defects are formed when variances are reduced by 
the continuing expansion of the Universe and a phase transition occurs. 



4.1 Non-thermal phase transitions 

Let us first describe ideas qualitatively. Let the infiaton oscillations which 
have amplitude 0(0) ~ Mp\ decay rapidly into some field X. In the instant 
process of decay the energy conserves and we have m^Mpj ~ k'^{X^), where 
on the left hand side of this equality we write the initial inflaton energy and 
on the right hand side we write the flnal energy stored in the X-fleld. The 
inflaton decays much faster than the system thermalizes, therefore typical 
momentum of X particles is of order of the inflaton mass, k ~ m^. We flnd 
(X^) ~ M|j. In thermal equilibrium with the temperature T we would have 
{X'^) = T^/12. In this respect the strength of fluctuations is the same as it 
would be in equilibrium with the Plankian temperature despite the fact that 
the real reheating temperature is much lower. Of course, this is an extreme 
estimate, mainly because the expansion of the Universe was neglected. The 
realistic {X"^) is shown in Fig. 6, but it can still exceed equilibrium temperature 
by many orders of magnitude. 

This effect is important for the behaviour of spontaneously broken symmetries 
at preheating. Let us consider the model in which in the vacuum state the 
symmetry is broken by an order parameter $. At the tree level this can be 
described by the potential 

= -i_i'^^^/2 + A$$V4- (21) 

The parameter is related to the symmetry breaking scale via /x^ = A$$q. If 
the flled $ and the product of inflaton decay, X, are coupled (with correspond- 
ing interaction term in the Lagrangian being aX^^^/2), at non-zero density 
of X-particles the effective mass of $ field changes to m^g = —fi'^ + a{X^). The 
symmetry is restored if the effective mass became positive, {X"^) > A$$q/q;. 

The real problem is complicated and model dependent. While some features 

can be anticipated and some quantities roughly estimated, the issue requires 
numerical studies. In recent papers [90-92] the defect formation and even the 
possibility of the first order phase transitions during preheating was demon- 
strated explicitly. 
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Fig. 7. Domain wall structure generated in simulation of preheating in the model 
Eq. (22) with M = 1, g = 0, Ref. [93]. 

4-2 Topological defects in simple models 

Let us consider for simplicity the system when one and the same field serves 
as the infiaton and the symmetry breaking parameter, $ = 0. We derive the 
set of models from the prototype potential 



The infiaton scalar field has M components, 0^ = YlfLi (ffi ^i^nd interacts 
with an A^-component scalar field X, = Y^f^iXf. For simplicity, the 
field X is taken massless and without self-interaction. The fields have minimal 
coupling to gravity in a FRW universe with a scale factor a{t). It is assumed 
that the infiaton oscillations start along 0i direction. In the effective mass of 
there are contributions from g'^{X^) as well as oc A ((50)^). 

We review several models, in the order of increasing complexity. Initial con- 
ditions in all cases correspond to the vacuum for fluctuations. System then 
evolved through particle creation, their rescattering, phase transitions and 
finally defects creation. 

4-2.1 Domain walls [93] 

Consider one component infiaton field, M = 1, and no X-fields. Relevant 
topological structures in the double-well potential are domain walls. Even 
in the model without fields X, fiuctuations of the field grow dramatically 




(22) 



19 




Fig. 8. String network generated in simulation of preheating in the model Eq. (22) 
with M = 2, g = (see Ref. [91]) is shown at two successive moments of time. 

(because of the self-interaction oc A, the zero mode can decay in the process 
20 —>■ 4:(j)] in addition to that the spinodal decomposition is effective in the 
present model). This leads to formation of a domain structure at certain values 
of A. The domain structure which emerged in such situation is shown in Fig. 7 



4.2.2 Strings [91] 

Again, consider the model without X-fields, but the inflaton has two compo- 
nents, M = 2. Relevant topological structures in the Mexican hat potential 
are strings. 

Fluctuations of the fields 0i and 02 during preheating restore symmetry along 
01 direction, but along 02 direction symmetry is broken. For some period of 
time the universe is divided into domains filled with the field 02 ~ ±v. Grad- 
ually the amplitude of fluctuations of the fields 0j decreases, and symmetry 
01 —01 also breaks down. At this moment the seed domain structure is 
transformed into a string network. Fig. 8 shows the string distribution in a 
simulation with the symmetry breaking scale v = 3 x 10^^ GeV, when a pair of 
"infinite" strings and one big loop were formed. Size of the box is comparable 
to the Hubble length at this time. 

The final result is the string formation, but the sequence of symmetry breaking 
patterns would be unusual for thermal phase transitions. For example, the 
0(2) symmetry was only partially restored. Another peculiarity of non-thermal 
phase transitions is the possible presence of oscillating zero mode during final 
stages of the phase transition. This makes probability of formation of long 
strings (as well as domains in the previous model) to be a non-monotonic 
function of f , at least for large v and in simplest models considered so far. 
In more complicated models inflaton oscillations can decay completely before 
the phase transition, and the phase transition can be even of the first order. 
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Defect formation in such models will be more robust and resembling more 
closely phase transitions in thermal background. 

4-2.3 First order phase transition [90] 

Here we consider one component inflaton field, M = 1 and several X-fields. 
In the case g^/X ^ 1 the phase transition can be of the first order. This 
can be expected by analogy with the usual thermal case [94]. The necessary 
conditions for this transition to occur and to be of the first order are as follows. 

(i) At the moment of phase transition a local minimum of the effective poten- 
tial should be at = 0, which gives g'^{X'^) > \v^. 

(ii) At the same time, the typical momentum of X particles should be 
smaller than gv. This is the condition of the existence of a potential barrier. 
Particles with momenta p < gv cannot penetrate the state with |0| ?a v, so 
they cannot change the shape of the effective potential at \(f)\ ^ v. Therefore, 
if both conditions (i) and (ii) arc satisfied, the effective potential has a local 
minimum at (j) — and two degenerate minima at ~ ±v. 

(iii) Before the minima at (f) ~ ±v become deeper than the minimum at (f) — 0, 
the inflaton's zero mode should decay significantly, so that it performs small 
oscillations near = 0. Then, after the minimum at |0| ^ v becomes deeper 
than the minimum at = 0, fluctuations of (p drive the system over the 
potential barrier, creating an expanding bubble. 

All these conditions can be met more easily at large g'^/X S> 1 and with several 
X-flelds. Results discussed below were obtained [90] on a 128^ lattice in the 
model with parameters g^/X = 200 and v = 0.7 x lO-^Mpi Ri 0.8 x lO^'' GeV, 
for a two-component X, with the expansion of the universe assumed to be 
radiation dominated (similar behaviour was observed in the model with single 
X field, i.e. = 1, as well). 

Time dependence of the zero mode 0o is shown in Fig. 9. Initially (po oscillates 
with a large amplitude 4> ^ v. If all fiuctuations were absent, the zero mode 
00, because of its dilution in the expanding universe, would at some time be 
unable to cross the potential barrier at v = and would start oscillations 
near one of its vacuum values, ±v. This would happen when the amplitude of 
the oscillations becomes smaller than \/2v. In Fig. 9 we see that the actual 
dynamics is completely different. The zero mode of the field continues to 
oscillate near = even when its amplitude becomes much smaller than 
V. In other words, the field oscillates on top of the local maximum of the 
bare potential. This can occur only because the effective potential acquires a 
minimum at = due to interaction of the field with {X"^). And amplitude 
of oscillations is decreasing much faster than it would decrease due to the 
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Fig. 9. Time dependence of the zero- momentum mode of (f) in units of its vacuum 
value for two runs with different realizations of random initial conditions for fluc- 
tuations in simulation of preheating in the model Eq. (22) with M = 1, ^ 7^ 0, Ref. 
[90]. 

simple expansion of the Universe. This happens because of continuing decay 
of the zero mode into X-particles. 

At still later times, r > 720 in Fig. 9, the zero mode of decays completely. 
This should be interpreted as restoration of the symmetry — > —0 by non- 
thermal fluctuations. Finally, at r > 860 (when X-fluctuations were diluted 
sufficiently by the expansion) a phase transition occurs and the symmetry 
breaks down. In runs with different realizations for initial ( "vacuum" ) ffuctu- 
ations the system ends up either in +v or in —v vacuum and the transition 
happens at different time moments. It was shown that transition is triggered 
in this model by a spontaneous nucleation of a bubble of the new phase and 
the bubble's subsequent expansion until it was occupying the whole integra- 
tion volume. The field configuration at the beginning of the phase transition 
is shown in Fig. 10, where the surface of the constant field = — 0.7f is plot- 
ted at the beginning of the phase transition. Inside the surface < —0.7^. 
(Minimum of the effective potential is shifted somewhat from the vacuum 
value at this time, see Fig. 9, because the contribution of ffuctuations is still 
significant) . 

Models that exhibit behaviour shown in Fig. 9 will lead to the domain structure 
surviving until present. This would be a cosmological disaster and such a class 
of models is ruled out [1]. 

However, different behaviour of the zero mode is observed at smaller values of 
/\ or (but still /\ S> 1). There, the phase transition can occur when the 
zero mode still oscillates near = with a relatively large amphtude. This 



22 



Fig. 10. Spontaneously nucleated bubble of the new phase at the initial stage of its 
expansion, Ref. [90]. 

is a new, specifically nonthermal, type of a phase transition. In such cases, 
bubbles of +f and -v phases will be nucleated in turn (most often at the 
maximum of amphtude of the zero mode when it is in the closest proximity 
to the top of the potential barrier), but their abundances need not be equal, 
and for certain values of the parameters one of the phases may happen not 
to form infinite domains. Such models are not ruled out and in fact may 
have interesting observable consequences, e.g. an enhanced (by bubble wall 
collisions) background of relic gravitational waves produced by the mechanism 
proposed in Ref. [18]. 

There are no reasons to doubt that in more complicated models non-thermal 
phase transitions may lead to creation of magnetic monopoles, or magnetic 
monopoles connected by strings (necklaces [43]), the latter being the most 
favourable topological defect candidate for explanation of UHECR [49] . 



5 Conclusions 

Next generation cosmic ray experiments, like the Pierre Auger Project [95], the 
High Resolution Fly's Eye [96], the Japanese Telescope Array Project [97], and 
the Orbiting Wide-angle Light-collector (OWL) [98], will tell us which model 
for UHECR may be correct and which has to be ruled out. 

Very weakly interacting superheavy X-particles with mx = (a few) ■ 10^^ GeV 
may naturally constitute a considerable fraction of Cold Dark Matter. These 
particles are produced in the early Universe from vacuum fluctuations dur- 
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ing or after inflation. Decays of X-particles may explain UHE cosmic rays 
phenomenon. Related density finctuations may have left an imprint in fluc- 
tuations of cosmic microwave background radiation if scalar X-particles with 
minimal coupling to gravity are approximately twice heavier than the inflaton 
and ~ 1- 

If UHE cosmic rays are indeed due to the decay of these superheavy particles, 
there has to be a new sharp cut-off in the cosmic ray spectrum at energies 
somewhat smaller than nix- Since the number density nx depends exponen- 
tially upon mx/mtj,, the position of this cut-off is fixed and can be predicted to 
be near 10^^ GeV, the very shape of the cosmic ray spectrum beyond the 

GZK cut-off being of quite generic form following from the QCD quark/gluon 
fragmentation. 

Very discriminating signature is related to anisotropy of cosmic rays. If par- 
ticles immune to CMBR are there, the UHECR events should point towards 
distant (i.e. beyond GZK sphere), extraordinary astrophysical sources [99]. If 
superheavy relic particles are in the game, the Galaxy halo will be reflected in 
anisotropy of the UHECR flux [100-102]. If neither wiU be true but UHECR 
will point instead to "local" galaxies with evidence for a central supermassive 
black hole, that would imply that the existence of a black hole dynamo is not a 
sufficient condition for the presence of pronounced jets and UHECRs are cre- 
ated by the remnants of dead quasars [32]. If none of the above will be true, 
arrival directions will be almost isotropic on large angular scales but UHE- 
CRs will cluster on small scales, then perhaps extragalactic magnetic fields 
are much stronger than previously thought [103] , or perhaps cosmic necklaces 
do exist. 

It is remarkable that we might be able to learn about the earliest stages of 
the Universe's evolution by studying UHECRs. Discovery of heavy relic X- 
particles will mean that the model of inflation is likely correct, or that at least 
standard early Friedmann evolution from the singularity is ruled out, since 
otherwise X-particles would have been inevitably overproduced [19,104]. 

We conclude that the observations of Ultra High Energy cosmic rays can probe 
the spectrum of elementary particles in the superheavy range and can give 
an unique opportunity for investigation of the earliest epoch of evolution of 
the Universe, starting with the amplification of vacuum finctuations during 
infiation through fine details of gravitational interaction and down to the 
physics of reheating. 



24 



6 Acknowledgements 



The work of V.K. was supported partially by the Russian Fund for Funda- 
mental Research grant 98-02-17493-a. 



References 



[I] Ya. B. Zeldovich, I. Yu. Kobzarev and L. B. Okun. Zh. Eksp. Teor. Fiz. 67 
(1974) 3 [Sov. Phys. JETP 40 (1974) 1 ]. 

[2] Ya. B. Zeldovich and M. Yu. Khlopov, Phys. Lett. 79B (1978) 239. 

[3] D. A. Kirzhnits, JETP Lett. 15 (1972) 529; D. A. Kirzhnits and A. D. Linde, 
Phys. Lett. 42B (1972) 471; Sov. Phys. JETP 40 (1974) 628; Ann. Phys. 101 
(1976) 195; S. Weinberg, Phys. Rev. D9 (1974) 3320; L. Dolan and R. Jackiw, 
Phys. Rev. D9 (1974) 3357. 

[4] Yu. Kobzarev, L. B. Okun, and M. B. Voloshin, Yad. Fiz. 20 (1974) 1229 [Sov. 
J. Nucl. Phys. 20 (1975) 644; S. Coleman, Phys. Rev. D15 (1977) 2922; S. 
Coleman and F. DeLuccia, Phys. Rev. D21 (1980) 3305; V. A. Berezin, V. A. 
Kuzmin, and I. I. Tkachev, Phys. Lett. 120B (1983) 91; J. Ipser and P. Sikivie, 
Phys. Rev. D30 (1984) 712; V. A. Berezin, V. A. Kuzmin, and I. I. Tkachev, 
Phys. Rev. D36 (1987) 2919; S. K. Blau, E. I. Guendelman, and A. H. Guth, 
Phys. Rev. D35 (1987) 1747. 

[5] A. A. Starobinskh, Phys. Lett. 91B (1980) 99. 

[6] A. H. Guth, Phys. Rev. D23 (1981) 347-356. 

[7] A. D. Linde, Phys. Lett. 108B (1982) 389; A. Albrecht and P. J. Steinhardt, 
Phys. Rev. Lett. 48 (1982) 1220. 

[8] A. D. Linde, Phys. Lett. 129B (1983) 177. 

[9] A. D. Linde, Particle Physics and Inflationary Cosmology (Harwood, Chur, 
Switzerland, 1990). 

[10] E. W. Kolb and M. S. Turner, The Early Universe, (Addison- Wesley, Reading, 
Ma., 1990). 

[II] V. F. Mukhanov and G. V. Chibisov, JETP Lett. 33 (1981) 532; S. W. 
Hawking, Phys. Lett. 115B (1982) 295 ; A. A. Starobinsky Phys. Lett. 117B 
(1982) 175; A. H. Guth and S.-Y. Pi, Phys. Rev. Lett. 49 (1982) 1110 . 

[12] G.F. Smoot et al., Ap. J. 396 (1992) LI; 

|http: / / www.gsfc.nasa.gov /astro / cobe / cobe_home.htmi . 

[13] |http:/ /map. gsfc.nasa.gov/ . 



25 



[14] |http:/ /astro. estec.esa.nl/Planck/ 



[15] |http:/ /www. sdss.org/ 



[16] L. Kofman, A. Linde, and A. A. Starobinsky, Phys. Rev. Lett. 76 (1996) 1011; 
I. I. Tkachev, Phys. Lett. B376 (1996) 35. 

[17] E. W. Kolb, A. D. Linde and A. Riotto, Phys. Rev. Lett. 77 (1996) 4290; 
E. W. Kolb, A. Riotto and I. I. Tkachev, Phys. Lett. B423 (1998) 348. 

[18] S. Y. Khlebnikov and I. I. Tkachev, Phys. Rev. D56 (1997) 653. 

[19] V. A. Kuzmin and I. I. Tkachev, JETP Lett. 68 (1998) 271. 

[20] K. Greisen, Phys. Rev. Lett. 16 (1966) 748; 

G. T. Zatsepin and V. A. Kuzmin, Pisma Zh. Eksp. Teor. Fiz. 4 (1966) 114. 

[21] N. N. Efimov et al., Proc. International Symposium on Astrophysicsl aspects 
of the most energetic cosmic rays, eds. M. Nagano and F. Takahara (1991, 
World Scientific, Singapore), p. 20; T. A. Egorov et. al., in: Proc. Tokyo 
Workshop on Techniques for the Study of Extremely High Energy Cosmic Rays, 
ed. M.Nagano (ICRR, U. of Tokyo, 1993); B. N. Afanasiev, Proc. International 
Symposium on Extremely High Energy Cosmic Rays: Astrophysics and future 
observations, ed. M. Nagano (Institute for Cosmic Rays Research, Tokyo, 
1996). 

[22] M. A. Lawrence, R. J. O. Reid, and A. A. Watson, J. Phys. G Nucl. Part. 
Phys. 17 (1991) 733; |http://ast.leeds . ac . uk / haver ah / hav- home . ht ml| . 

[23] D. J. Bird et. al, Phys. Rev. Lett. 71 (1993) 3401; Astroph. J. 424 (1994) 491; 
441 (1995) 144. 



[24] N. Hayashida et. al, Phys. Rev. Lett. 73 (1994) 3491; |http: / /icrsun.icrr.u- 
I t okyo . ac . j p / as / pr o j ect / agasa. ht mj| . 

[25] V. S. Berezinsky, S. V. Bulanov, V. A. Dogiel, V. L Ginzburg, and V. S. 
Ptuskin, "Astrophysics of cosmoc Rays", Elsevier, 1990. 

[26] P. Bhattacharjee and G. Sigl, |astro-ph /98110li ; V. Berezinsky, [astro- 



ph/9811268 



[27] Takeda M. et al., Phys. Rev. Lett. 81 (1998) 1163. 

[28] M. A. Ruderman and Sutherland, Ap. J. 196 (1975) 51; R. D. Blandford and R. 
Znajek, MNRAS 179 (1977) 433; A. M. Hillas, Ann. Rev. Astron. Astrophys. 
22 (1984) 425. 

[29] C. A. Norman, D. B. Melrose, and A. Achterberg, Ap. J. 454 (1995) 60. 
[30] J. W. Elbert and P. Sommers, Ap. J. 441 (1995) 151. 



[31] A. Dar, |astro-ph/9901005 



[32] E. Boldt and P. Ghosh, |astro-ph/9902342 



26 



[33] S. Coleman and S. L. Glashow, [hep-ph/OSOSlie . 

[34] G. R. Farrar, Phys. Rev. Lett. 76 (1996) 4111; D. J. Chung, G. R. Farrar, and 
E. W. Kolb, Phys. Rev. D57 (1998) 4696. 

[35] S. Raby, Phys. Lett. B422 (1998) 158. 

[36] M. B. Voloshin and L. B. Okun, Sov. J. Nucl. Phys. 43 (1986) 495. 

[37] T. J. Weile r, |hep-ph/971043l| ; D. Fargion, B. Mele, and A. Salis, |asfao^ 
I ph/9710029| ; S. Yoshida, G. Sigl, and S. Lee, Phys. Rev. Lett. 58 (1998) 55. 

[38] G. Gelmini and A. Kusenko, |hep-ph/990235^ . 

[39] J. Bahcall and E. Waxman, |hep-ph/9902383 . 

[40] W. H. Kinney and A. Riotto, [hep-ph/9903459| ; J. Adams and S. Sarkar, talk 
presented at the ICTP Workshop on physics of relic neutrinos, Trieste, Italy, 
1998. 

[41] C. T. Hill, D. N. Schramm and T. P. Walker, Phys. Rev. D36 (1987) 1007. 

[42] J. H. MacGibbon and R. H. Brandenberger, Nucl. Phys. B331 (1990) 153; P. 
Bhattacharjee and N. C. Rana, Phys. Lett. B246 (1990) 356. 

[43] V. Berezinsky, X. Martin and A. Vilenkin, Phys. Rev. D56 (1997) 2024; V. 
Berezinsky and A. Vilenkin Phys. Rev. Lett. 79 (1997) 5202. 

[44] R. L. Davis and E. P. S. Shellard, Phys. Lett. B209 (1988) 485; R. 
Brandenberger et al, Phys. Rev. D54 (1996) 6059; C.J.A.P. Martins and E. 
P. S. Shellard, |hep-ph/9806480| . 

[45] C. T. Hih, Nucl. Phys. B224 (1983) 469; P. Bhattacharjee and G. Sigl, Phys. 
Rev. D51 (1995) 4079. 

[46] T. W. Kephart and T. J. Weiler, Astropart. Phys. 4 (1996) 271. 

[47] V. A. Kuzmin and V. A. Rubakov, Talk at the Workshop "Beyond the Desert" , 
Castle Rindberg 1997; Phys. Atom. Nucl. 61 (1998) 1028 [Yad. Fis. 61 (1998) 
1122]. 

[48] V. Berezinsky, M. Kachehiess and A. Vilenkin, Phys. Rev. Lett. 79 (1997) 
4302; 

[49] V. Berezinsky, P. Blasi, and A. Vilenkin, Phys. Rev. D58 (1998) 103515. 
[50] M. Birkel and S. Sarkar, Astropart. Phys. 9 (1998) 297. 

[51] K. Benakh, J. Elhs, and D. V. Nanopoulos, Phys. Rev. D 59 (1999) 047301 ; K. 
Benakh, Phys. Lett. B447 (1999) 51; D. V. Nanopoulos, |liep-ph/9809"54e . 

[52] T. Han, T. Yanagida and R.J. Zhang, Phys. Rev. D58 (1998) 095011. 

[53] K. Hamaguchi, Y. Nomura, and T. Yanagida, Phys. Rev. D58 (1998) 103503; 
Phys. Rev. D59 (1999) 063507. 



27 



[54] M. Kachelriess, |hep-ph/9806322| . R. Dick, N. Eschrich and M. Gaul, |ep| 
ph/9808245|; G. K. Leontaris and J. Rizos, |hep-ph/990T098|; K. Hamaguchi, 



K. I. Izawa and Y. Nomura, |hep-ph/9903207 . 

[55] L. Kofman, A. Linde, and A. Starobinsky, Phys. Rev. Lett. 73 (1994) 3195; 

[56] S. Y. Khlebnikov and I. I. Tkachev, Phys. Lett. B390 (1997) 80. 

[57] S. Y. Khlebnikov and I. 1. Tkachev, Phys. Rev. Lett. 79 (1997) 1607. 

[58] T. Prokopec and T. G. Roos, Phys. Rev. D55 (1997) 3768; B. R. Greene, T. 
Prokopec, and T. G. Roos, Phys. Rev. D 56 (1997) 6484. 

[59] D. J. H. Chung, E. W. Kolb and A. Riotto, Phys. Rev. D59 (1999) 023501. 

[60] V. A. Kuzmin and I. I. Tkachev, |hep-ph/9809547 . 

[61] J. Ellis, J. Kim and D. V. Nanopoulos, Phys. Lett. B145 (1984) 181. 

[62] E. W. Kolb, D. J. H. Chung, and A. Riotto, Phys. Rev. Lett. 81 (1998) 4048; 
|hep-ph/9809455 . 

[63] S. G. Mamaev, V. M. Mostepanenko, and A. A. Starobinskii, ZhETF 70 (1976) 
1577 [ Sov. Phys. JETP 43 (1976) 823]. 

[64] S. G. Mamaev, V. M. Mostepanenko, and V. M. Frolov, Yad. Fiz. 23 (1996) 
1118, ibid 26 (1997) 215. 

[65] L. Parker, Phys. Rev. 183 (1969) 1057; A. A. Grib and S. G. Mamaev, Yad. 
Fiz. 10 (1969) 1276 [Sov. J. Nucl. Phys. 10 (1970) 722 ]; Ya. B. Zel'dovich and 
A. A. Starobinslii, Zh. Exp. Theor. Fiz. 61 (1971) 2161 [Sov. Phys. JETP 34 
(1972) 1159]; L. P. Grishchuk, Sov. Phys. JETP 40 (1975) 409. 

[66] A. A. Grib, S. G. Mamaev, and V. M. Mostepanenko, Quantum Effects in 
Strong External Fields, [in Russian], Atomic Energy Press, Moscow (1980). 

[67] N. N. Bogolubov, Sov. Phys. JETP 7 (1958) 51. 

[68] N. A. Chernikov and E. A. Tagirov, Ann. Inst. Henri Poinare 9 A, 109 (1968); 
E. A. Tagirov, Ann. Phys. 76 (1973) 561; T. S. Bunch and P. C. W. Davies, 
Proc. R. Soc. A360 (1978) 117. 

[69] A. D. Linde, Phys. Lett. 116B (1982) 335; A. A. Starobinsky, Phys. Lett. 117B 
(1982) 175; A. Vilenkin and L. H. Ford, Phys. Rev. D26 (1982) 1231; B. Allen, 
Phys. Rev. D32 (1985) 3136; C. M. Massac and and C. Schmid, Ann. Phy s. 
231 (1994) 363; H. F. Muller and C. Schmid, |gr-qc/941202l| ; |gr-qc/941202^ . 

[70] A. A. Starobinsky and J. Yokoyama, Phys. Rev. D50 (1994) 6357. 



[71] D. Lyth and D. Roberts, piep-ph/9609441 . 

[72] K. Enqvist, K. W. Ng and K. A. Olive, Nucl. Phys. B303 (1988) 713. 



28 



[73] A. A. Starobinsky and V. Sahni, in: Proc. 6th Soviet Gravitation Conf., ed. V. 
N. Ponomarev (MGPI Press, Moscow, 1984); G. Efstathiou and J. R. Bond, 
Mon. Not. R. Astron. Soc. 218 (1986) 103; H. Kodama and M. Sasaki, Int. J. 
Mod. Phys. A2 (1987) 491; A. A. Starobinsky, Sov. Astron. Lett. 14 (1988) 
166. 

[74] A. Linde and V. Mukhanov, Phys. Rev. D56 (1997) 535. 
[75] P. J. E. Peebles, |astro-ph/9805T9^ ; |astro-ph/98052T^ . 

[76] L. H. Ford, Phys. Rev. D35 (1987) 2955; M. Joyce an d T. Prokopec, Ph ys. 
Rev. D57 (1998) 6022; P. J. E. Peebles and A. Vilenkin, |astro-ph/9810509| ; G. 
Felder, L. Kofman, and A. Linde, hep-ph /9903350| . 

[77] J. Preskill, M. B. Wise, and F. Wilczek, Phys. Lett. 120B (1983) 127; L. F. 
Abbott and P. Sikivie, Phys. Lett. 120B (1983) 133. 

[78] J. Ellis, D. V. Nanopoulos and M. Quiros, Phys. Lett. B174 (1986) 176. 

[79] I. I. Tkachev, Sov. Astron. Lett. 12 (1986) 305; I. I. Tkachev, Phys. Lett. B191 
(1987) 41; I. I. Tkachev, Phys. Lett. B261 (1991) 289. 

[80] J. Garcia-Belhdo and A. Kusenko, |hep-ph/9812"53^ . 

[81] T. W. Kephart and T. J. Weiler, Phys. Rev. D52 (1995) 3226. 

[82] C. J. Hogan and M. J. Rees, Phys. Lett. B205 (1988) 228; E. Kolb and I. I. 
Tkachev, Phys. Rev. Lett. 71 (1993) 3051; Phys. Rev. D 49 (1994) 5040; Phys. 
Rev. D50 (1994) 769; Ap. J. 460 (1996) L25. 

[83] A. D. Dolgov and A. D. Linde, Phys. Lett. B116 (1982) 329; L. F. Abbott, 
E. Farhi, and M. B. Wise, Phys. Lett. B117 (1982) 29; A. Albrecht, P. J. 
Steinhardt, M. S. Turner, and F. Wilczek, Phys. Rev. Lett. 48 (1982) 1437; A. 
D. Dolgov and D. P. Kirilova, Sov. Nucl. Phys. 51 (1990) 172; 

[84] J. Traschen and R. Brandenberger, Phys. Rev. D 42 (1990) 2491; Y. Shtanov, J. 
Traschen, and R. Brandenberger, Phys. Rev. D 51 (1995) 5438; D. Boyanovsky, 
H. J. de Vega, R. Holman, D. S. Lee, and A. Singh, Phys. Rev. D 51 (1995) 
4419; M. Yoshimura, Prog. Theor. Phys. 94 (1995) 8873; D. Kaiser, Phys. Rev. 
D 53 (1996) 1776. 

[85] L. Kofman, A. Linde, and A. A. Starobinsky, Phys. Rev. D 56 (1997) 3258. 
[86] S. Y. Khlebnikov and I. I. Tkachev, Phys. Rev. Lett. 77 (1996) 219. 
[87] D. Polarski and A. A. Starobinsky, Class. Quant. Grav. 13 (1996) 377. 
[88] D. T. Son, Phys. Rev. D54 (1996) 3745. 

[89] D.J.H. Chung, |hep-ph/9809'489t G. Felder, L. A. Kofman, and A. D. Linde, 
|hep-ph/981228g| . 

[90] S. Khlebnikov, L. Kofman, A. Linde, and I. Tkachev, Phys. Rev. Lett. 81 
(1998) 2012. 



29 



[91] I. Tkachev, S. Khlebnikov, L. Kofman, and A. Linde, Phys. Lett. B440 (1998) 
262. 



[92] S. Kasuya and M. Kawasaki, Phys. Rev. D58 (1998) 083516. 

[93] S. Khlebnikov, L. Kofman, A. Linde, and L Tkachev, unpublished; I. Tkachev, 
talk at 18th Texas Symposium on Relativistic Astrophysics, Chicago, IL, 15-20 



Dec 1996, |hep-ph/9701376 



[94] D.A. Kirzhnits and A.D. Linde, Ann. Phys. (N.Y.) 101 (1976) 195. 

[95] J. W. Cronin, Nucl. Phys. B (Proc. Suppl.) 28 (1992) 213; The Pierre Auger 
Observatory Design Report (2nd edition), March 1997; M. Boratav, for Pierre 
Auger Collaboration, The Pierre Auger Observatory Project: An Overview., 
Proc. of 25th International Cosmic Ray Conference, Durban, v. 5, p. 205 



(1997); http://www.auger.org 



[96] S. C. Corbato et al, Nucl. Phys. B (Proc. Suppl.) 28 (1992) 36; 

|http:/ /bragg. physics. adelaide.edu.au/astrophysics/HiRes. html 



[97] M. Teshima et al, Nucl. Phys. B (Proc. Suppl.) 28 (1992) 169; 
|http://www-ta.icrr.u-tokyo.ac.jp| . 

[98] J. F. Ormes et al., Proc. 25th International Cosmic Ray Conference (Durban, 
1977), eds: M. S. Potgieter et al.. Vol. 5, 273; 

|http: / /Iheawww. gsfc.nasa.gov /docs/ gamcosray /hecr / OWL / . 



[99] G. R. Farrar and P. L. Biermann, Phys. Rev. Lett. 81 (1998) 3579. 
[100] S. L. Dubovskii and P. G. Tinyakov JETP Lett. 68 (1998) 107. 
[101] V. Berezinskii and A. A. Mikhailov, |astro-ph/98102771 . 



[102] M. Takeda et al, |astro-ph/990223£ . 



[103] M. Lemoine, G. Sigl, and P. Biermann, |astro-ph/9903124| ; D. Ryu and P. L 



Biermann, Astron. Astrophys. 335 (1998) 19. 
[104] I. I. Tkachev, talk at the conference on Particle Physics and the Early Universe 



(COSMO 98), Monterey, CA, 15-20 Nov 1998; |hep-ph/9901428 



30 



